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1 Introduction 

The amount of isoscalar tensor states claimed to exist experimentally [jl| has reached a point 
where naive interpretation of these states becomes perilous. This is particularly distressing 
in light of the fact that the tensor glueball, a degree of freedom beyond conventional mesons, 
is the second lowest glueball predicted by lattice QCD [§, ||, ^ and should be manifested 
in this multiplicity of states. There is a need to bring some order by sorting out well 
understood states. What is sorely needed is a model-independent theoretical tool. 

The need for model-independence of the theory is especially prevalent in light of the fact 
that the recently claimed 02 (1660) |^, |, ^, |] has a mass which confounds traditional, 
and often reliable, potential models of radially excited P-wave mesons. For example, the 
difference between the first radially excited and ground state isovector J^^ = 2++ (tensor) 
states were predicted to be 510 MeV in a relativized quark model 0, while the experimental 
value is ~ 340 MeV. 

In this paper we present a mass-matrix analysis of considerable, although not total, gen- 
erality. Schwinger-type mass formulae [^ are derived. It is assumed that the mesons 
and glueball mix only via meson-glueball coupling, with no direct meson-meson coupling. 
At any stage of the analysis we restrict ourselves to a finite number of mesons (with one 
glueball) and mixing with hypothetical four-quark states is not taken into account. 

The term isoscalar "mesons" shall refer to the partners of the light quark isovector mesons, 
each of which has an ss partner. The isovector mesons will be given labels like P-wave, F- 
wave or hybrid meson, indicating the dominant component in a quark model interpretation 
of the state. However, the mass matrix analysis does not assume the P-wave, F-wave or 
hybrid meson nature of a state, nor that it is an unmixed quark model state. 

2 Masses 

The tensor sector has a few salient features which simplify the analysis. The first excited 
lattice QCD tensor glueball is 1.85 ± 0.20 times the mass of the scalar glueball P, and its 
effect on the experimental spectrum can hence safely be neglected. This means that we 
can restrict consideration to the low-lying mesons and one primitive (bare) glueball. The 
mass of the glueball is reliably estimated by using M(2++)/M(0++) = 1.39 ± 0.04 or 



1.42 ±0.06 0, in combination with the average lattice QCD value M(0++) ^ 1.6 GeV ^, 
to be around 2.2 GeV 0]. It is instructive to obtain the lower limit on the tensor glueball 
mass allowed by lattice QCD. With 1.5 GeV the lower limit for the scalar glueball mass, 
one obtains a tensor glueball mass ~ 1.35 M(0^+) ~ 2.0 GeV. This limit will be employed 
later on. 

The isovector tensor mesons should act as beacons for the mass scales of various nonets. 
Unfortunately, only the 02(1320), which we will take to fix the mass of the primitive nn 
ground state P-wave state, IP, is well established [jl|. There is recent evidence for 02(1660) 
at 1660±40 MeV or 1660±15 MeV §, for an 02(1600-1700) §, evidence at ARGUS for 
a mainly 2++ state at 1.7 GeV § and for an 02 at 1752 ± 21 ± 4 MeV [|. The 02(1660) is 
taken to fix the mass of the primitive nn first radially excited P-wave state, 2P. Additional 
evidence for the presence of a 2P nonet is provided by its 1++ partners. The Oi(1700) was 
claimed by BNL [12| and a similar signal was seen at VES @, 13|. Recently, weak evidence 



for /i(~ 1700) was reported |T^. 

There is also some recent evidence for isovector tensor states at 2060 ± 20 MeV and 1990^30 



MeV [0, signalling the 3P and IF nonets. Except for these isovector tensor states, the 
reasons for expecting the 3P and IF nonets in a similar mass region are as follows. There 
are /o(2010), /o(2060) |I|], and /j(2100), with J most likely 0, at 2115 ± 15 ± 15 MeV §. 
There are recent indications of an oq at 2025 ± 30 MeV ||T5[. These J = states signal 
P-wave mesons, since neither the ground state nor the excited scalar glueball is expected 
in this mass region ||^. An oi at 2100 ± 20 MeV also indicates 3P. Given the experimental 
mass splitting between the 2P and IP nonets noted earlier one also expects the 3P level in 
this mass region. We shall take the primitive nn 3P level to be at 2.05 GeV. The IF nonet 
is signposted by the O4(2050), /4(2050), i\:4(2045), the recently reported 03 at 1860 ± 20 



MeV H or 2070 ± 20 MeV p[ and /a at 2000 ± 40 MeV |T1 or 1950 ± 15 MeV p 
There are recent indications from VES ^ that the mass of the O4(2050) is 1944 ± 8 ± 50 
or 1950 ± 20 MeV. We place the primitive nn IF state at 1.94 GeV and the primitive ss 
state higher by twice the difference between the fC4(2045) and VES' O4(2050) masses, i.e. 
at 2.15 GeV. Variation of these masses is discussed in Appendix 0. 



There is some evidence for an isovector tensor state at 2265 ± 20 MeV |T5[, signalling the 



2F or 4P nonets. The presence of both nonets is indicated by an /i at 2340 ± 40 MeV [|14 



an oi at 2340 ± 40 MeV ^, and an /o at 2335 ± 25 MeV ^, which can be 4P but not 



2F; or /4(2300), 7^3(2320) [0], an a^ at 2300 ± 20 MeV ig, an /a at 2280 ± 30 MeV 0, 
and as at 2310 ± 40 MeV, which can be 2F but not 4P [|. 

It is clear that there is no evidence for overpopulation of levels for isovector tensors, implying 
that there is no need to introduce a hybrid meson level0 up to ~ 2.3 GeV. For isodoublet 
tensors, there is in fact an underpopulation: only the well-established ii'2(1430) from the 
IP nonet, and the marginal -ft'lllQSO) are known ||l|. 

To the contrary, there are 13 isoscalar tensor mesons up to ~ 2.3 GeV listed by the Particle 
Data Group, with 6 well-established^ '^. One expects a glueball, and the IP, 2P, 3P and 
IF nonets in this mass region, yielding 9 states, and possibly nn 4P and 2F in addition, 
giving 11 states. There is hence an overpopulation of experimental isoscalar tensors, albeit 
not for the well-established ones. 

Since the IP, 2P, 3P and IF nonets are expected below ~ 2.3 GeV, our analysis can safely 
be restricted to a 9 x 9 mass matrix. There is the possibility of the 4P and 2F mesons 
contaminating results at the upper end of our simulation, at ~ 2.3 GeV, which is also 
investigated. 



3 5x5 mass matrices 

The mixing of a glueball and n pairs of isoscalar mesons is described by the following mass 

matrix, motivated in Appendix |B|, which is diagonalized by the masses of (2n+ 1) physical 

^/4(2300) and 04 may be members of the IH nonet, although the nonet appears to be more high-lying, 

as signalled by the 06(2450) and /6(2510) |]. 

^ The possibility of a tensor hybrid meson in the mass range up to ^ 2.3 GeV cannot be excluded 

theoretically. Beyond the early MIT bag model estimates, constituent gluon models have estimated a 

tensor hybrid mass, most recently at 1.6 — 1.8 GeV [Eol. Lattice QCD splittings of hybrid levels indicate 

that at least for bb hybrids, the tensor hybrid is degenerate with the lightest hybrids within errors [[41[ . 

However, adiabatic lattice QCD and flux-tube models do not find tensors on the lowest hybrid adiabatic 

surface. Also, tensor mesons are associated with 0+"*" hybrids in bag, constituent gluon and flux-tube 

models and adiabatic lattice QCD. There is no indication of an overabundance of isovector scalar states. 

^Taking both /j(1710) and /j(2220) to have 2++ components. 



states: 



G gi giV2 §2 g2\f^ ■ ■ ■ 9n 9nV^ \ 

gi Si ■ ■ ■ 

giV2 Ni ■ ■ ■ 

^2 52 ■ ■ • 

^2^2 A^2 • • • 



(1) 



gn ■ ■ ■ Sn 
\ gnV2 ■ ■ ■ Nn J 
=^ diag {hi, h2, hs ..., /i2n, h2n+i)- 

G and 5", A^ stand for the mass of the primitive glueball, and ss and nn = {uu + dd)/ \/2 
mesons, respectively, the subscript indicating the number of the nonet the state belongs 
to. hi stand for the masses of the physical states, g^ are the glueball-meson couplings that 
have dimensionality (mass), in accord with the dimensionality of the diagonal entries of 
(H). In what follows, we restrict ourselves to the case where the quantities in ([^) are real 
numbers^. 

Applying the techniques of ref. |T^, one can obtain 2n pairs of relations for the coupling 
in terms of the primitive and physical masses: 
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gi 
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(2) 



Each pair of these relations represents a Schwinger-type mass formula. Hence, for (2n + 
1) X (2n + 1) mass matrix (||) one has n Schwinger mass relations. These n formulae, 
together with the trace condition for the mass matrix (|l]). 



G + S1 + N1 + S2 + N2 



Sn + Nn = hi + h2 + h3 + . . . + h2n + /i2n+l, 



(3) 



constitute n + 1 mass relations for the mixing of a glueball and n meson nonets. It is 
clear that solving such a system of n + 1 mass relations can lead to unphysical solutions. 



*gi and —gi gives the same eigenvalues, so we always choose gi non-negative. 



e.g., solutions that correspond to all or some of the couplings being imaginary numbers. 
Obviously, such solutions will not correspond to the initial mass matrix (|l]), and hence 
should be rejected. 

It is difficult to find all the solutions of the Schwinger equations for a 9 x 9 mass matrix 
of the form (|I]) numerically. Hence we take the approach to solve the Schwinger equations 
for a 5 X 5 mass sub-matrix, which involves only the primitive and physical masses, and 
then reconstruct the couplings. As all 5 x 5 sub-matrices are found, we then obtain the 
corresponding 9x9 mass matrix. 

In this analysis, we pursue the following strategy: 

(i) We start with the 5x5 sub-matrix for the glueball and 2P and IF nonets, by fixing 
the masses of primitive nn for 2P, and nn and ss for IF, and obtain the primitive 2P ss 
mass. We then consider another 5x5 sub-matrix for the glueball and 2P and 3P nonets, 
with fixed: both nn and ss masses of the 2P nonet from the previous simulation, and the 
primitive nn mass for the 3P nonet. In both 5x5 sub-matrix simulations we obtain all 
the solutions of the Schwinger equations (H). 

(ii) Input: For the two 5x5 sub-matrix simulations, we fix the following values of the 
primitive masses (in GeV): A^ = 1.66 for 2P, A^ = 1.94, 5 = 2.15 for IF and A^ = 2.05 for 
3P. We also take one of the physical states to have a mass in agreement with one of the 
glueball candidates (which we review in the next section), and the other three physical states 
to have masses in agreement with three states among /2(1565), /2(1640), /j(1710), /2(1810), 
/2(1950), /2(2010), /2(2150), /j(2220), /2(2300) and /2(2340), excluding the state already 
chosen for the physical glueball. Output: We then solve the system of three equations 
(two Schwinger formulae and the trace condition) for three unknowns: G, S for 2P and 
the remaining fifth physical mass for the first 5x5 simulation, and G, S for 3P and the 
remaining fifth physical mass for the second 5x5 simulation. We require that the fifth 
physical mass from each 5x5 simulation is among the physical states mentioned above. 

(iii) Since we take the /2(1275) and /2(1525) as the established ground state IP tensor 
mesons, we incorporate them later in the full 9x9 mass matrix analysis. 

(iv) We discard the possibility that /2(1420) exists. Although claimed by a number of old 
experiments in a variety of production processes, recent experiments do not confirm its 
existence. This is most vividly illustrated by its observation in (mostly) double Pomeron 



exchange in pp -^ pfln^ii )ps at ^/s = 63 GeV |^, |T6]. Recent examination of the same 
reaction does not see any evidence for /2(1420) |T^. 

We admit the following criteria for holding physical solutions and separating out non- 
physical ones: 

(i) The output fifth physical mass lies within a mass range allowed by data for one of the 
experimental candidates. 

(ii) The mass of the primitive glueball satisfies G >2 GeV. 

(iii) In all the cases when a primitive ss mass is to be obtained, it is higher than the 
corresponding nn mass, and the ss — nn mass splitting is consistent with the quark model 



motivated estimate 200 ± 50 MeV 11 



4 9x9 mass matrices 

Various physical states have been suggested as tensor glueball candidates in the literature: 

/j(2220): The /j(2220) is strongly produced in J/ip radiative decay, and not seen in 77 
collisions, suggesting glueball character if J = 2 |^, |I8|. The flavour democratic decay 
pattern and small total width of /j(2220) is also cited as evidence for its glueball nature 

m. 



/2(2150): This was suggested in ref 

The nearness of the mass of /j(2220) and /2(2150) to the tensor glueball mass predicted 
by lattice QCD is often cited as evidence for their glueball nature 0. 

/2(1950): The pr dependence of the pp central production of /2(1950) is consistent with its 
glueball character according the Close-Kirk glueball fllter [^ ^. However, it was admitted 
that the structure seen in central production may represent more than one resonance pT|. 



/j(1710): The glueball nature of this state is suggested by its px dependence in central 
production ^l| and its production in "glue-rich" pp annihilation [|l8l, although its strong 



production in J /%p radiative decay is consistent with expectations for qq if J=2 [|T^ . 

For each 5x5 case, we take one of the above four glueball candidates to be one of the flve 
physical states. 



Having completed the double 5x5 sub-matrix analysis and fixing A^ = 1.318 for IP ||T|, the 
full 9x9 mass matrix is now recovered by solving the Schwinger equations @ exactly, using 
the solutions obtained for the two 5x5 matrices as initial values for the search routine. 
We therefore do not obtain all of the 9x9 mass matrix solutions, but only the ones similar 
to the ones found formerly with the two 5x5 matrices. 

There are two different solutions, which are almost identical with respect to the physical 
masses. Particularly, the primitive glueball masses are consistent with 2.0 — 2.1 GeV pre- 
dicted by recent models |^, ^ and the lattice QCD predictions mentioned earlier. The 
couplings are in the range 30 — 120 MeV for the various nonets. These values are similar to 



43 ± 31 MeV predicted in lattice QCD for ground state isoscalar scalars [^. We find that 
the physical masses are insensitive to changes in the input, but that the valence content is 
more sensitive: especially for states at similar masses to where the parameters are changed, 
and for small valence components (see Appendix]^). 

(i) For the first solution, the /j(2220) turns out to be the physical glueball. For this 
solution, the initial 9x9 mass matrix is (shown are the values of the primitive masses and 
couplings rounded to the second decimal digit; all values are given in GeV) 

^ 2.10 0.03 0.03^2 0.04 0.04^2 0.09 0.09^2 0.12 0.12^2 \ 

0.03 2.28 

0.03^2 2.05 

0.04 2.15 

0.04^2 1.94 , (4) 

0.09 1.84 

0.09^2 1.66 

0.12 1.55 

\ 0.12v^ 1.318 J 

the physical masses are 



2.29, 2.23, 2.14, 2.04, 1.93, 1.82, 1.64, 1.52, 1.2J 



(5) 



and the valence content of the physical states is 
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(6) 



(ii) For the second solution, the the physical glueball is distributed, with /2(1810) containing 
the largest component. Although the highest mass state appears to have the largest glueball 
component, we shall see in section ^ that the content changes as more high mass states are 
introduced. For this solution, the initial 9x9 mass matrix is 
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(7) 



2.38, 2.23, 2.12, 2.01, 1.95, 1.82, 1.63, 1.52, 1.2J 



(8) 



and the valence content of the physical states is 
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(9) 



5 Discussion of mass matrix results 



Both solutions have the following similarities: 

(i) The physical states are /2(1270), /;(1525), /2(1640), /2(1810), /2(1950), /2(2010), /2(2150), 
/j(2220) and either the /2(2300)//2(2340) or /2(2340) as solutions in the first and second 
cases, respectively]^. We never found the /2(1565) and /j(1710). The reason why /2(1810) 
is found instead of these resonances is because the primitive 2P ss is required to 250 ± 50 
MeV from the input 02(1660) mass. 

(ii) The valence content has almost entirely the same signs between the various components, 
the only exception being different signs for the two dominant meson components in /2(1950). 

(iii) The physical mesons have a substantial glueball content, contrary to naive expectations, 
with the exception of /2(1950) in solution 2. It has been argued phenomenologically that 
experimental data demand physical mesons with appreciable glueball content [|^. This 
would, for example, explain why /2(2010), /2(2300) and /2(2340) were observed in the OZI 
forbidden process vrp — ^ 0077,, and would suggest that several tensor mesons should be 



produced in glue-rich processes. For example, /2(1270) was observed in gluon fusion [^ 



The small glueball component in /2(1950) in solution 2 is apparently in contradiction with 
the Close-Kirk filter. 



'''Within experimental mass uncertainty [El, we cannot distinguish between /2(1525) or /2(1565). In 

TTTT is consistent with experiment if it is taken to have the valence 



section 7.1, we show that /2(1525) 

content of the second state, which we thus identify as /2(1525). 
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(iv) The /2(1270), /2(1525) and /2(1640) are composed of more than 90% of the expected 
primitive state. 

(v) For the IP nonet, ^ + A^ = 1.55 + 1.318 = 2.868 GeV is consistent with 2M{K^) = 
2.858 ±0.01 GeV H. 

The sohitions differ as follows: 

(i) In the first solution all physical mesons have one component which has a valence content 
of larger than 90%, i.e. the state is dominantly a specific primitive state. For this solution, 
valence components of physical mesons greater than 10% only occur within two primitive 
states of the dominant primitive state. The second solution does not behave in this way. 

(ii) For the first solution, the physical glueball has substantial valence content in all the 
meson states it couples to, contrary to the second solution. 

(iii) In the first solution the couplings decrease with increased radial excitation (from IP to 
3P), as one would naively expect [^; while for the second solution the couplings remain 
approximately the same, as expected from Regge theory [^. An advantage of our method 
is that all information on couplings are predictions, using (0), once the masses are known. 

(iv) The first solution has the slight disadvantage that the primitive ss-nn mass splitting 
~ 180 MeV in the 2P nonet is a little small, versus an agreeable ~ 290 MeV for the second 
solution. 

(v) For the second solution, it is seen in (|^) that both the /2(1950) and /2(2150) have large 
nn and ss components, and neither of them are in destructive interference. Hence, both 
of these states will have many different decay modes, which is in excellent agreement with 
data (these different decay modes are amongst the main reasons for each of these states 
to be chosen as the tensor glueball candidate by different groups). Although this feature 
is absent for these states for the first solution (§)-(§]), where /2(1950) is mostly nn and 
/2(2150) mostly ss, many of the observed decay modes arise from connected decay of both 
uu and ss components, so that this avenue to distinguish between solutions may not be 
definitive. 
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5.1 /2(1565) and /j(1710) 

It is possible that /2(1565) and /2(1640) are aspects of the same state, which would remove 
one extra state. We have nevertheless attempted to find solutions where /2(1565) and 
/j(1710) are physical states in addition to the states (H),®, by adding another nonet to 
form an 11 x 11 matrix. Of course, one can insert primitive states at these masses with 
an unrealistically small glueball-meson coupling and an unrealistic ss-uu mass splitting of 
150 MeV and obtain a solution. However, no realistic solutions are found. 

We have neglected the mixing of mesons with decay channels thoughout, since it is believed 



to produce only tiny mass shifts [^. However, scenarios where there is large coupling to 
decay channels, e.g. ujoj, pp, K*K* and 00, have been advanced by various authors. 

/2(1565) decays to pp and ujuj and has an abnormally small branching ratio to tttt and rji] 
|l|]. This, together with the nearness of /2(1565) to the pp and ujuj thresholds has lead to 
suggestions that /2(1565) is a pp molecule or a baryonium state |2^. Also, /2(1565) may 
be the isoscalar partner of the isotensor tensor enhancement X(1600) [|l|, which, if it is 
resonant, must be a degree of freedom beyond glueballs and (hybrid) mesons. 

A more modest suggestion is that the mass of the 2P nn state found in our formalism is 
shifted downward by the pp and ujuj thresholds, which the ^Pq model predicts it to couple 
strongly to pO|. 



/j(1710) has been suggested as a K*K* molecule |^. However, it is not well established 



that a J = 2 component exists. BES separated both J = and J = 2 components, with 



the tensor state having mass 1697 MeV and a width of 176 MeV [^. However, recent 
evidence supports only the J = component |jl^ ^ . 



6 13 X 13 mass matrix 

Once the 9x9 mass matrix is fixed, one can easily add extra meson nonets to it. We add 
the primitive nn and ss masses of the 2F nonet at 2.3 and 2.5 GeV, respectively, to the 
mass matrices @,(0). Similarly, the primitive states of the 4P nonet are added at 2.35 and 
2.55 GeV, to yield a 13 x 13 mass matrix. The physical states in Eqs. d^),® are required 
to be among the physical states, i.e. both /2(2300) and /2(2340). 
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The result is a 13 x 13 "counterpart" to each 9x9 matrix (|^),(0), called solutions la and 
2a. The primitive states in common have the same couplings and primitive masses, and 
similar valence content (with the same signs). The valence content of a given primitive 
state tends to decrease from the 9x9 counterpart, since the physical state is spread over 
more primitive states. 

Remarkably, the ratios of valence contents of the 13 x 13 solutions and their 9x9 coun- 
terparts remain extremely similar (~ 1%), except for the components in the 9x9 matrix 
which has similar mass to the new components being added. This means that for low-lying 
states, there is usually no need extend the number of primitive components in order to 
study decay. 

We also find two new solutions, called lb and 2b, since they are respectively similar to the 
9x9 solutions 1 and 2. They have, however, no 9 x 9 counterparts. Solution 2b is displayed 
in Appendix |^. 

In all 13 X 13 solutions, three new, experimentally undiscovered, physical states appear at 
masses beyond the /2(2340). The dominant glueball component in solutions 2a and 2b is 
found in one of the three new high mass states. 

We note that because the number of nonets stable under decay is expected to be finite due 
to pair creation in QCD [Q, the largest mass matrix that need to be analysed is finite. 

7 Decays 

In order to calculate the decay of a physical state to an exclusive final state it is necessary 
to add the decay amplitudes of all its primitive components, weighted by their valence 
content. This is demonstated in Appendix 0. 

The decay amplitudes of the primitive components will be calculated in the ^Pq model, 
meaning that pair creation is with vacuum quantum numbers and decays proceed via a 
connected quark diagram. Unless otherwise noted, the decays are calculated using the 
"relativistic" phase space convention and parameters of ref. |^. Another convention is 
"mock meson" phase space with the parameters of ref. |3^ . In the mass matrix analysis, 
specific quark model identifications were not assumed for the various components. However, 
to calculate the decays, the quark model content indicated by the label of a component will 
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be assumed. 

The preceding mechanism whereby the physical glueball decays via primitive meson compo- 
nents, which is closely related to the primitive glueball decaying via intermediate primitive 
mesons [^ , provides the first theoretical understanding of the scalar glueball decay pattern 
found in lattice QCD ^. 

When one combs the experimental data available on isoscalar tensor mesons more massive 
than the /2(1525), there is very little robust quantitative information available that can 
directly be compared to theory. The most restrictive datum appears to be the ratio of 
widths of /j(2220) to tttt and KK, which we shall analyse below. Of the qualitative data 
available, the observation or non-observation of various isoscalar tensor mesons in 00 is 
of particular interest, as the 00 decay can arise for connected decay only from the ss 
components of states. We shall not make exhaustive decay predictions, but restrict to the 
cases mentioned in a genuine attempt to confront our picture with experiment. However, 
we first analyse the OZI forbidden decay /2(1525) -^ tttt which is zero in models where the 
state has only one valence component. 

7.1 /2(1525) -^ TTTT 

For the 13 X 13 solutions la, lb, 2a and 2b we find r(/2(1525) ^ tttt) = 1.6(1.4), 1.2(1.1), 1.0 
(0.9), 0.7(0.7) MeV, with relativistic phase space^ listed first. Solution 2b is in best agree- 
ment with the experimental value 0.60 ± 0.12 MeV. The two 9x9 solutions gives the same 
results as their 13 x 13 counterparts. As seen in (|D,(||), the valence content of the /2(1525) 
is such that its IP and 2P nn components are in destructive interference (they have opposite 
signs), which will result in the suppression of the tttt decay mode of this state. Furthermore, 
for components higher in mass than 2P, the valence content has the same sign as the 2P 
component, leading to further suppression. To be specific, we illustrate this for solution 2b. 
r(/2(1525) -^ tttt) = 14 and 4 MeV if only IP, and IP and 2P components are included. 
The width remains above 1.5 MeV as long as not all of the IP, 2P, IF and 3P components 
are included. 



^ In accordance with ref. |Q we use a slightly higher pair creation constant for low-mass states. For 
r(/2(1275) —* tttt), with /2(1275) purely IP nn, this gives 160 MeV, in perfect agreement the experimental 



157 ±4 MeV [§. 
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We have thus provided the first quantitative understanding of the process /2(1525) -^ tttt. 
This demonstrates that the techniques of both the mass matrix and ^Pq decay analysis yield 
predictions consistent with experiment, motivating their continued use. Remarkably, the 
decay /2(1525) -^ vrvr can only be understood when at least four different components of 
/2(1525) are included. It is also apparent that there is no need to postulate a non-connected 
decay mechanism, whereby primitive ss components would directly decay to tttt. Because 
/2(1525) is dominantly ss, such processes must be small indeed. 

7.2 R = r(/j(2220) -^ 7r+7r-)/r(/j(2220) ^ K+K ) (J = 2) 

The main distinguishing characteristic of the /j(2220) is its remarkably narrow total width 
of 23^7 MeV ^ for a state that can decay via numerous decay modes [^. Not even the 
existence of /2(2220) is well-established 0] as the narrow peak sit on a variety of empirical 
backgrounds, for which there is no explanation ||3^ , so that the peak might be a statistical 
fiuctuation. Moreover, broader states in the same mass region have been reported: JETSET 
sees an /a at 2231 ± 2 MeV with a width of 70 ± 10 MeV p8l . An /a at 2240 ± 40 MeV with 



a width of 170 ± 50 MeV H, and at 2210 ± 45 MeV with width 260 ± 45 p] have also 
been reported. It is possible to understand current data on /j(2220) if one does not take 
it to be narrow |^. If the /j(2220) is not narrow, none of the "indicators" of its glueball 



nature can be sustained, for example its non-observation in 77 collisions [^ simply follows 
from its wideness, and its coupling to gluons in J/\E' radiative decay becomes compatible 
to conventional mesons if J = 2 IH 



As seen in (^),(^, the largest non-glue components of the /j(2220), 3P and IF ss, are in 
destructive interference, and remain so in the 13 x 13 solutions. This also tends to be true 
for the largest nn components in the 13 x 13 solutions. One may think this will result in 
the suppression of the decay modes of this state, making its total width consistent with the 
tiny experimental value. Evaluating the total width of /j(2220), with J = 2, to tttt and 
KK, for all the 13 x 13 solutions, and phase space conventions, we obtain 20 — 150 MeV. 
It is evident that the tiny total width can not be sustained in our model, and it likely is 
to be a challenge to any model in which the physical glueball has non-negligible mixing 
with mesons [^. The individual partial widths to vrvr and KK are also inconsistent with 



experimental bounds that assume a narrow fj (2220) 
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It is often claimed that /j(2220) has a flavour democratic decay pattern expected for a 



pure glueball |T9|, whereby R = 1 without phase space included, and R = 1.7 with phase 
space included. However, naive flavour factors give that a pure nn and ss state should 
have, without phase space included, R = 4,0 respectively. Thus a mixture between nn and 
ss can also look flavour democratic. For the 13 x 13 solutions la, lb, 2a and 2b we find 
R = 0.6, 0.7, 0.4, 0.5 respectively, independent of phase space conventions. The two 9x9 
solutions gives the same results as their 13 x 13 counterparts. Although these values of 



R do not represent flavour democratic decay, they are all consistent with experiment |]19 
which possesses large error bars. This is true for solutions la and lb where /j(2220) is the 
physical glueball, and for the other solutions. 

When decays are to final S-wave mesons, i.e. 7r,ri,K, p,uj,K*,ri or </>, one almost always 
finds that the decay amplitudes decrease sharply as the decaying component is progressively 
radially excited . The same is true as the decaying component is orbitally excited. This has 
the consequence that although a physical state may have a dominantly excited component, 
its decay dominantly proceeds through a lower excited component. This means that naive 



quark model calculations that assign a single component to an excited state pO| , p5| might be 
completely unreliable. One would a priori expect this situation to be worst in J^'-' sectors 
where low-lying glueballs are present, i.e. J^*" = 0'^"*",2++ and 0~+ 0. We illustrate the 
phenomenon by analysing /j(2220) -^ tttt and KK for the 13 x 13 solutions. Although 
f J (2220) is never dominantly IP, this contribution is always one of the dominant ones. 
One tends to finds that half of the width can be found by including only the IP and 2P 
contributions, even though the state may be dominantly IF and 3P. 

7.3 Decays to 00 

We shall analyse the decay of various resonances to (jxf), in an attempt to understand the 
data, which claim that /2(2010), /2(2300) and /2(2340) have been seen in (jxp in irp collisions 
||1|]. There is also preliminary evidence for an /2 at ~ 2231 ± 2 MeV in (fxp [^. We bear in 
mind that the production process can alter conclusions made based on studying decays. 

The results are in Table [^. From phase space considerations, it is especially surprising 
that it is possible for /2(2150) to have a smaller width to (jxp than /2(2010). This is even 
more surprising, given that /2(2010) is dominantly nn and /2(2150) is dominantly IF ss. 
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13 X 13 solution 


/2(2010) 


/2(2150) 


/j(2220) 


/2(2300) 


/2(2340) 


la 


3(3) 


7(5) 


17(12) 


23(16) 


7(4) 


lb 


3(2) 


5(4) 


12(9) 


17(12) 


8(5) 


2a 


18(15) 


3(3) 


0.6(0.5) 


4(3) 


9(6) 


2b 


12(11) 


1.7(1.4) 


0.7(0.5) 


7(5) 


2(2) 



Table 1: Decay widths to 00 in MeV for various solutions. Relativistic phase space is 
used first, with mock meson phase space in brackets. The two 9x9 solutions give almost 
identical results to their 13 x 13 counterparts. /2(2010) is taken to have a mass at the upper 
end of the experimental range [Q . 

The resolution of the paradox is that IF and 2F does not decay to 00 in the dominant 
S-wave in the ^Pq model 0]. We note that solutions 2a and 2b are consistent with states 
unambiguously observed in 00. 

8 Salient features 



We showed that Schwinger-type mass formulae can be obtained when we restrict to glueball- 
(hybrid) meson mixing. With some physical isovector and isoscalar masses known, these 
formulae can predict unknown masses and couplings. The utility of this new analysis tech- 
nique was demonstrated in the tensor sector. 

It has been shown that in order to understand the decay /2(1525) -^ tttt, one has to consider 
more than the IP nn component. This implies that the use of a 2 x 2 mixing formula where 
the physical states are linear combinations of nn and ss components can be inadequate. 

In our approach the physical glueball is a priori narrower than mesons due to the large 
glueball component, which is taken not to decay. However, as shown for the /j(2220), 
the physical glueball is not unusually narrow, because of the presence of significant meson 
components, contrary to the perturbative QCD claim that glueball mixing with IF mesons 
is tiny [^. Also, there is no reason to expect a flavour democratic decay pattern for the 
physical glueball. Experimentally, this is already clear for the scalar glueball 0, |11|. It is a 
common myth that primitive glueballs are narrow and decay flavour democratically. These 
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notions arise from perturbative QCD, which, paradoxically, has been argued to be valid for 



the tensor glueball but not for the scalar glueball |39|. However, glueballs beyond 3 GeV 



can be narrow [^. A more reliable glueball signature may be glue-rich production. 
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A Appendix: 13 x 13 mass matrix 

The 13 X 13 mass matrix for solution 2b is 



i 2.25 


0.08 


0.08 V 


^ 0.08 


0.08 %/2 


0.13 


0.13V2 


0.14 


0.14V 


^ 0.08 


0.08 V2 


0.12 


0.12V2 \ 


0.08 


2.55 



































0.08 ^/2 





2.35 
































0.08 








2.5 





























0.08 V2 











2.3 


























0.13 














2.27 























0.13\/2 

















2.05 




















0.14 




















2.15 

















0.14\/2 























1.94 














0.08 


























1.95 











0.08^2 





























1.66 








0.12 
































1.55 





i^ O.12V2 



































1.318 J 

(A.l) 


rhe physical masses 


are 




















2.67, 


2.53, 


2.47, 


2.34, 


2.29, 2.23 


, 2.12 


, 2.01, 


1.95, 


1.81, 


1.64, 1.52, 1.28 




(A.2) 
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and the valence content of the physical states is 



0.65 


0.45 


0.23 


0.32 


0.20 


0.21 


0.19 


0.18 


0.18 


0.07 


0.07 


0.07 


0.08 \ 


0.19 


-0.84 


0.12 


0.46 


0.09 


0.09 


0.07 


0.07 


0.06 


0.02 


0.02 


0.02 


0.03 


0.27 


-0.29 


0.25 


-0.82 


0.18 


0.18 


0.12 


0.12 


0.10 


0.04 


0.04 


0.04 


0.04 


0.11 


-0.04 


-0.88 


-0.05 


0.37 


0.23 


0.07 


0.09 


0.06 


0.02 


0.02 


0.02 


0.02 


0.08 


-0.02 


-0.15 


-0.03 


-0.79 


0.57 


0.06 


0.08 


0.05 


0.02 


0.01 


0.01 


0.01 


0.23 


-0.06 


-0.21 


-0.07 


-0.36 


-0.69 


0.24 


0.43 


0.16 


0.07 


0.05 


0.04 


0.04 


-0.17 


0.03 


0.08 


0.04 


0.11 


0.15 


-0.44 


0.83 


-0.19 


-0.08 


-0.04 


-0.04 


-0.04 


-0.18 


0.03 


0.06 


0.03 


0.07 


0.09 


0.74 


0.17 


-0.54 


-0.26 


-0.06 


-0.05 


-0.04 


-0.02 


0.00 


0.00 


0.00 


0.00 


0.01 


0.03 


0.01 


-0.38 


0.92 


-0.01 


-0.00 


-0.00 


-0.43 


0.05 


0.09 


0.05 


0.10 


0.12 


0.32 


0.17 


0.64 


0.24 


-0.33 


-0.20 


-0.15 


-0.20 


0.02 


0.03 


0.02 


0.03 


0.04 


0.09 


0.05 


0.13 


0.05 


0.92 


-0.28 


-0.10 


-0.22 


0.02 


0.03 


0.02 


0.03 


0.04 


0.08 


0.05 


0.10 


0.04 


0.18 


0.93 


-0.18 


-0.23 


0.01 


0.02 


0.01 


0.03 


0.03 


0.05 


0.04 


0.07 


0.03 


0.07 


0.10 


0.96 1 



(A.3) 



Solution 2b is similar to the 9x9 solution 2 in Eqs. (0)"®' ^^^ ^^ ^°^ ^^^ 13 x 13 
counterpart of it, as evidenced by the different couplings of the common primitive states 
and the different primitive glueball mass. 

We now study the stability of the solution under parameter changes. When we change 
the primitive IF nfi and ss masses upwards by 50 and 40 MeV respectively, reflecting 
our lack of knowledge of these parameters, the physical masses all remain consistent with 
experiment, except for /j(2220), which is 6 MeV higher than the experimental mean [Q. 
Valence contents of states far away in mass from IF remain essentially constant. The largest 
changes are found for /2(2010), where valence contents less than 0.2 change on average by 
70%. The dominant content changes by 1% and the next most dominant by 20%. A further 
change of the 2P ss mass downwards by 30 MeV yields the largest changes for the valence 
content of /2(1950), by similar amounts as before. 

B Appendix: Motivation for mass matrix 

If an appropriate basis is chosen, the hamiltonian for n isovector states up to a certain mass 
can be taken to be a diagonal n ^ n matrix M = diag(A'^i, A^2, • • • , Nn)-, where the entries 
are real and positive. These entries are identified with the masses of the physical states in 
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the experimental isovector spectrum. Note that no assumption is made about the nature 
of the states, e.g. conventional or hybrid meson or four-quark state. In a world of only 
u, d quarks, this work assumes degeneracy of isovector and isoscalar states, motivated in 
ref. [jlO|- The hamiltonian for the isoscalar states is simply the (n + 1) x (n + 1) matrix 

G V2g , , ^ 

B.l 

V2g^ N ^ 

where a new state, which cannot exist in the isovector sector, the glueball, has been added. 
G is the primitive glueball mass, and g is a. n-dimensional row of (real) couplings of either 
uu or dd. These couplings are 0{^=) in the large number of colours N^ expansion of QCD. 
Of course, any number of extra glueballs can in principle be added ||10||. 



If the strange quark is also incorporated, the isoscalar matrix becomes (|lD. Here the matrix 
element between nn states i and j {i ^ j) is zero because of the diagonality of Af in ( [B.l|) . 
The strange quark states are taken to be heavier analogues of the light ones, so that the 
mixing between ss states i and j {i ^ j) is zero, and the diagonal entries contain the masses 
Si of the strange quark states. The remaining possible mixing is between nn state i and 
ss state j, which is 0{^) [^ and was found tiny in recent lattice QCD simulations p3 



and is neglected in this work. The strange quark states are assumed to couple in the same 
way to the glueball as uu and dd, which is the SU(3) limit. For ground state isoscalar 
scalar mixing, lattice QCD obtains the ratio of uu and dd to ss glueball coupling to be 



1.198 ± 0.072 H, while the SU(3) limit is 1. 



C Appendix: Decay formalism 

The mass matrix in Eq. |1] can be viewed as forming part of a hamiltonian H that describes 
an effective theory. This part of the hamiltonian is diagonalized to yield the physical states. 
One can then a posteriori add to the hamiltonian a part that describes coupling to a decay 
channel. 



H = g*Gg + J2is:S^s, + n:Nin,) + iJ2g,g*{s, + V2n,) + c.c.) + {Yg* + J2iY^s:+7>mBC) 

i=l i=l 1=1 

(C.l) 
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where G, Ni and Si are the masses and gi the couphngs in Eq. |1]; g,ni,Si the corresponding 
primitive glueball, nn and ss meson fields; and ■j^y'Ji'y'yi represent the couphngs of the 
respective primitive states to the decay channel field (BC). Spin indices for the primitive 
tensor states have been suppressed in H, e.g. g*Gg stands for g*^^Gg'^^, where g^'^ is a 
symmetric and traceless Lorentz tensor. 

Now write the hamiltonian in shortened notation 



H = m^Mui + m^^{BC) m = (c/, Si,ni . . . s„, n„)^ -f = {-ff^ , Yi, li ■ ■ ■ In^nf 

(C.2) 
where A1 is defined as the matrix (|l|). Since A4 is real and symmetric it can be diagonalized 
by use of the (real) orthogonal {Q^^ = Q'^) valence content matrix Q, to yield the diagonal 
mass matrix of (real) eigenvalues (physical masses) A4o = QAiQ^^. The eigenvectors 
(physical states) are m = Qm, where m denotes the primitive states. The j^^ row in 
Q gives the valence content of the j*^ physical state. Several practical examples of the 
diagonalization procedure can be found in this work. For example, for the Al in (^), JAd 
is the matrix with diagonal entries (§), and Q is (^). 

The hamiltonian becomes 

H = m^MDm + m^^-f{BC) (C.3) 

The first term was discussed in detail in ref. [^. The second term shows clearly that in 
order to calculate the decay amplitude of a physical state to (BC), it is necessary to add 
the decay amplitudes of all its primitive components, weighted by their valence content. 



This was used, but not explicitly demonstrated in ref. pT]] . 
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